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I.  INTRODUCTION 


In  a  recent  study,  120- mm  gun  firings  were  performed  to  establish  the  ballistic 
contributions  due  to  the  combustible  cases  used1.  Measured  velocities  and  pressures  were 
compared  with  matching  interior  ballistic  code  simulations  with  the  hope  of  quantifying  those 
contributions.  Previously,  for  firings  for  which  the  charge- mass-tc>  projectile-mass  (C/M)  ratio 
was  low,  IBHVG2,  a  lumped-parameter  interior  ballistic  code,  has  provided  good 
comparisons1’2.  Also  in  the  past,  the  NOVA  family  of  two-phase  interior  ballistic  codes,  of 
which  XNOVAKTC  (XKTC)  is  the  latest  version,  has  been  shown  to  be  able  to  simulate  firings 
with  very  good  suocess1’2’3’4.  In  this  recent  study,  however,  for  which  the  C/M  was  about  unity, 
there  was  a  wide  disparity  between  the  predictions  of  XKTC  and  those  of  IBHV02  for  the  same 
nominal  data  base. 

In  the  previous  work1,  it  was  shown  that  XKTC  could  be  made  to  mimic  the 
experimental  120-mm  gun  firings  quite  well.  For  IBHVG2,  however,  that  was  not  the  situation. 
For  one  case,  IBHVG2  gave  a  calculated  maximum  breech  pressure  that  was  42  MPa  higher 
than  that  predicted  by  XKTC  for  the  same  nominal  data  base.  It  was  found  that  a  difference  of 
14  MPa  could  be  attributed  to  the  fact  that  IBHVG2  does  not  model  the  projectile  boattail 
intrusion,  and  that  3  MPa  each  could  be  attributed  to  the  fact  that  IBHVG2  neither  models 
flamespreading  nor  intergranular  stress.  The  major  difference  in  the  predicted  maximum 
pressures,  however,  was  attributed  to  the  simple  physics  used  in  the  derivation  of  the  pressure 
gradients  allowed  by  IBHVG2.  In  XKTC,  the  axial  pressure  gradient  is  calculated  from  first 
principles  and  analytic  correlations;  in  IBHVG2,  only  analytical  pressure  gradient  relations  due 
to  Lagrange  and  due  to  Pidduck-Kent5  are  available. 

In  this  report,  several  ballistic  parameters  which  might  affect  the  pressure  gradient, 
especially  for  large  C/M  ratios,  are  examined  in  some  detail.  The  objective  of  this  study  was  to 
determine  the  physics  that  must  be  included  in  the  analytic  gradient  equation,  so  that  the 
predictions  of  lumped-parameter  codes  can  be  improved. 

II.  INITIAL  COMPARISONS  BETWEEN  EXPERIMENT  AND  XKTC 

The  120-mm  experimental  gun  firings  are  described  in  a  previous  paper1.  The  gun  tube 
was  instrumented  with  five  pressure  gages  in  the  chamber:  two  at  95  mm,  one  at  286  mm 
(midchamber)  and  two  at  489  mm  from  the  rear  face  of  the  tube.  There  were  also  seven 
downbore  gages,  situated  at  768  mm,  1048  mm,  1530  mm,  2292  mm,  3054  mm,  3816  mm,  and 
4578  mm  from  the  rear  face  of  the  tube.  Assuming  that  the  gages  at  95  mm  determine  the 
breech  pressure  accurately,  that  the  base  of  the  projectile  before  it  moves  is  located  at  541  mm 
from  the  rear  face  of  the  tube,  and  that  the  pressure  measured  as  each  downbore  gage  is 
uncovered  is  the  projectile  base  pressure  at  that  time,  one  can  use  the  data  contained  in 
Reference  1  to  calculate  the  ratio  of  the  breech  pressure  to  the  projectile  base  pressure  for 
several  discrete  values  of  projectile  travel.  Table  1  presents  the  average  pressure  ratios  for  the 
three  gun  firings  which  were  performed  with  no  cases  and  average  pressure  ratios  for  the  three 
gun  firings  which  were  performed  with  inert  cases. 
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Tkble  1.  Experimental  Ratios  of  Breech  Pressure  to  Projectile 
Base  Pressure  at  Several  Discrete  Values  of  Projectile  Travel 
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The  data  points  for  the  ratios  of  the  breech  to  projectile  base  pressure  versus  travel  from 
Tfcble  1  are  plotted  against  calculated  curves  of  these  ratios  from  XKTC  in  Figure  1.  Both  the 
caseless  gun  firing  series  and  the  inert  case  gun  firing  series  are  shown.  Breech  and  base 
pressure  curves  are  added  to  assist  in  interpreting  the  ratio  data.  The  XKTC  calculations  were 
performed  using  measured  values  wherever  possible  and  reasonable  values  for  all  other  input 
data.  These  simulations  gave  excellent  agreement  with  measured  pressure-time  curves  and 
pressure  difference  curves.  The  agreement  between  the  measured  values  and  the  calculated 
values  for  pressure  ratios  for  all  but  the  last  two  points  is  very  good.  The  reason  for  the  last 
points  not  fitting  well  is  not  known.  In  any  case,  this  good  agreement,  at  least  for  the  major 
portion  of  the  ballistic  event,  should  make  XKTC  a  useful  tool  for  studying  the  gradient 
phenomenology.  With  this  understanding,  in  this  report,  XKTC  calculations  have  been  assumed 
to  furnish  the  "correct*  answers  with  which  to  compare  lumped-parameter  calculations. 
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Figure  1.  Comparisons  Between  Measured  and  XKTC-Calculated  Ratios  of  Breech  Pressure  to 
Base  Pressure. 


III.  FIRST  MODEL  COMPARISONS 

As  we  have  shown  that  that  XKTC  could  be  used  for  accurate  simulations  of  gun  firings, 
we  sought  to  perform  calculations  with  XKTC  specifically  simplified  so  that  IBHVG2 
calculations  could  be  compared  to  them.  To  that  end,  the  first  series  of  calculations  performed 
with  XKTC  utilized  data  bases  with  (a)  7-perforated  granular  propellant  evenly  distributed  along 
the  length  of  the  chamber,  (b)  the  chamber  diameter  equal  to  the  bore  diameter,  (c)  zero  barrel 
resistance,  (d)  a  flat-based  projectile,  (e)  nominal  heat  loss  to  the  chamber  walls,  and  (f)  a 
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propellant  ignition  temperature  equal  to.  the  ambient  temperature^  ao  that  the  prapeilam  Was  all 
ignited  at  the  sun  of  calculations.  A  typical  XKTC  data  base  used  in  this  study  is  included  as 
Appendix  A.  The  companion  IBHV02  daU  base  is  included  u  Appendh  B. 

»  '  • 

*  .  fV  ,*  «  |  »  ,• 

Guns  with  two  different  chamber  volumea  were  simulated,  one  with  a  chambar  Volume 
of  96312  enr  and  one  with  a  chamber  tribune  of  96322  cm r.  TWo  different  guns  were 
simulated  in  order  to  determine  whether  the  obeervad  effects  were  a  faction  of  the  physical 
size  of  the  weapon.  The  other  weapon  parameters  that  were  associated  with  sech  of  thiaae  two 
guns  are  are  drawn  in  Thble  2.  For  each  of  the  two  weapon*  two  different  projectile  Weights 
were  used  to  produce  the  two  different  C/M  ratios  for  which  calculations  were  performed. 


Table  2.  Weapon  Parameters  for  the  First  Series  of  Calculations 


Chamber  volume 

98323  cm3 

98.322  cm3 

Travel 

4.572  m 

1.88  m 

Propellant  mass 

9.8  kg 

0.098  kg 

Projectile  masses 

9.8  kg  and  393  kg 

0.098  kg  and  0392  kg 

Bore  diameter 

127.0  mm 

28.65  mm 

The  effects  -of  propellant  bum  rate  and  the  maximum  chamber  pressure  were  also 
examined.  Table  3  shows  the  parameters  which  were  varied: 

Table  3.  Other  Parameters  for  the  First  Series  of  Calculations 


Burning  Rate  (p  in  MPa) 

1.10519  p1 0  mm/s 

0.408451  p®  8  mm/s 

Max  Pressure 

172  MPa 

345  MPa  517  MPa 

The  propellant  thermochemistry  for  all  calculations  is  shown  in  Table  4. 

% 

Table  4.  Propellant  Thermochemistry 


Impetus 

1136  J/g 

Co  volume 

0.976  cm3/g 

Gamma 

1.23 

Flame  1'emperature 

3141  K 

Molecular  Weight 

23.0  g/g-mde 

Density 

1.66  g/cm3 
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Table  5.  Comparimm  of  Calculated  Maximum  Breech  Ptbmutm 
md  Munte  Vdoddea  for  XKTC  and  IBHV02 

XKTC  EBHVQ2 
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The  comparison  calculations  performed  with  XKTC  and  with  IBHVG2  (using  the 
Lagrange  gradient  model  and  the  Pidduck-Kent  gradient  model)  are  summarized  in  Table  5  for 
the  variations  described  above.  In  the  table,  the  "Ch  Vol"  is  the  chamber  volume  of  the 
particular  weapon,  and  “BR"  is  the  propellant  burning  rate,  where  ”1.0"  implies  1.10519  p10 
mm/s  and  *0.8"  denotes  0.408451  p0  8  mm/s.  The  maximum  breech  pressure  and  muzzle  '/elocity 
with  their  associated  times  are  given  in  the  table.  For  each  horizontal  line  on  this  table,  the 
propellant  web  in  the  XKTC  calculation  was  varied  until  the  desired  peak  pressure  was  achieved; 
the  propellant  length  was  maintained  between  2  and  3  times  the  outer  propellant  diameter. 
Then,  with  the  same  final  propellant  dimensions,  the  associated  IBHVG2  Lagrange  and  the 
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IBHVG2  Pidduck-Kent  calculations  were  performed 

Close  inspection  of  this  table  reveals  that  IBHVG2  agrees  with  XKTC  very  closely  when 
C/M  is  025,  for  both  the  Lagrange  gradient  and  the  Pidduck-Kent  gradient,  for  both  weapons 
and  for  both  propellant  burr  rates.  However,  when  C/M  is  1.0,  agreement  is  not  as  good 

Tins  first  comparison  series  of  calculations  with  XKTC  and  IBHVG2  investigated  the 
influence  of  C/M,  chamber  volume,  propellant  burning  rate  and  maximum  breech  pressure  on 
ballistic  performance.  Figures  2  through  5  have  been  generated  from  XKTC  calculations  for  the 
weapon  with  the  large  chamber  volume.  Figure  2  shows  plots  of  the  mean  to  projectile  base 
pressures  for  C/M  equal  to  1.0,  for  increasing  maximum  breech  pressures,  and  for  different 
burning  rates.  Figure  3  shows  plots  of  the  breech  to  mean  pressure  ratios  for  C/M  equal  to  1.0, 
for  increasing  maximum  breech  pressures,  and  for  different  burning  rates.  Figure  4  shows  plots 
of  the  mean  to  projectile  base  pressures  for  C/M  equal  to  0.25,  for  increasing  maximum  breech 
pressures,  and  for  different  burning  rates.  Figure  5  shows  plots  of  the  breech  to  mean  pressure 
ratios  for  C/M  equal  to  0.25,  for  increasing  maximum  breech  pressures,  and  for  different 
burning  rates.  In  each  of  these  figures,  tlie  plots  on  the  let!  have  increasingly  higher  maximum 
breech  pressures  and  were  performed  with  a  burning  rate  of  1.10519  p10  mm/s,  while  the  plots 
on  the  right  have  the  same  increasingly  higher  maximum  breech  pressures  and  were  performed 
with  a  burning  rate  of  0.408451p°  8  mm/a  The  shape  and  magnitude  of  the  plots  for  the  smaller 
chamber  volume  were  the  same  for  the  corresponding  C/M,  burning  rates  and  pressures  except 
the  time  scale  and  travel  scale  were  about  1/3  that  of  the  larger  chamber  volume  calculations, 
so  they  have  not  been  included  in  this  report.  As  an  aid  in  interpreting  the  ratios,  the  baJe  and 
breech  pressures  are  plotted  on  each  graph. 
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Figure  5.  XKTC-Calculated  Breech/Mean  Pressure  Ratio  Curves  for  a  C/M  of  0.25 

The  ratio  plots  in  Figures  2  through  5  illustrate  the  complex  nature  of  the  relationship 
between  the  breech,  mean  and  base  pressures.  The  general  nature  of  the  ratio  curves  is  to  be 
constant  near  1.0  for  some  length  of  time,  followed  by  a  first  rise,  a  drop-off,  a  second  rise  and 
in  some  case};  a  second  drop-off  and  another  rise.  The  second  drop-off  seems  to  be  caused  by 


the  slivering  phase  of  the  propellant  burning,  as  it  occurs  just  as  XKTC  indicates  that  slivering  is 
taking  place.  Further,  for  low  pressure  calculations,  for  which  the  propellant  does  not  sliver 
before  projectile  exit,  the  second  drop-off  does  not  occur.  Calculations  done  with  a  single 
perforation  grain  also  did  not  show  a  second  drop-off.  The  first  rise  and  drop-off  may  be  caused 
by  the  rarefaction  wave  caused  by  the  projectile  motion,  the  first  peak  corresponding  to  the  time 
that  the  rarefaction  wave  reaches  the  breech  face  and  the  first  minimum  corresponding  to  the 
time  the  reflected  wave  reaches  the  base  of  the  projectile.  In  all  of  these  XKTC  calculations, 
we  observed  that  an  ullage  region  opens  up  between  the  projectile  base  and  the  front  end  of  the 
propellant  bed.  In  this  ullage  region,  the  pressure  drop  per  unit  of  distance  across  this  ullage 
region  is  much  smaller  than  across  the  propellant  bed.  Gough5  has  speculated  that  the 
formation  of  an  ullage  region,  as  observed  in  the  XKTC  calculations,  and  the  discontinuity  in  gas 
velocity  at  the  propellant/ullage  boundary  may  contribute  to  the  undulatory  shape  of  the  mean 
to  base  pressure  ratio. 

There  remained  the  question,  however,  of  what  physical  processes  were  responsible  for 
the  greater  than  10%  maximum  breech  pressure  differences  observed  between  real-world 
simulations  with  XKTC  and  the  equivalent  IBHVG2  calculations,  the  differences  that  motivated 
this  study. 


IV.  SECOND  MODEL  COMPARISONS 

The  second  set  of  model  comparisons  included  calculations  at  both  high  and  low  C/M, 
because  C/M  was  implicated  in  the  lack  of  agreement  between  XKTC  a  id  IBHVG2  in  the  first 
comparisons.  The  Lagrange  gradient  was  used  for  all  IBHVG2  calc  jlations.  Several  other 
parameters  were  also  varied  to  determine  their  contribution  to  the  problem.  For  these 
calculations,  we  modeled  the  weapon  with  the  9832.2  cm3  chamber  volume,  and  we  used  a 
burning  rate  exponent  of  1.0.  We  varied  the  barrel  resistance  for  both  XKTC  and  IBHVG2 
calculations,  and  for  XKTC,  introduced  flamespreading  and  added  chambrage  (the  necking  down 
of  the  chamber  to  the  bore  diameter). 

Since  a  constant  hamber  volume  was  desired,  chamber  length  had  to  change  when 
chambrage  was  introduced.  The  chambrage  was  modeled  as  a  truncated  cone  whose  length  was 
76.2  mm,  which  required  the  chamber  length  be  reduced  from  776.22  mm  to  541.02  mm,  with 
the  radius  of  the  breech  end  of  the  chamber  up  to  the  beginning  of  the  chambrage  being  76.91 
mm.  Baseline  and  chambrage  configurations  are  illustrated  in  Figure  6. 

Flamespreading  in  XKTC  is  convectively  driven,  with  the  initial  stimulus  provided  by 
some  level  of  modeling  of  igniter  functioning.  In  these  calculation,  the  igniter  was  described  as 
venting  over  the  rear  304.8  mm  of  propellant  bed,  causing  the  entire  propellant  bed  to  be  ignited 
within  2  ms.  The  resistance  was  modeled  with  a  linearly  interpolated  table  of  travel  versus 
resistive  pressures.  The  resistance  started  at  0.690  MPa  at  0-mm  travel,  remained  constant  for 
6.35  mm,  rose  to  6.90  MPa  at  12.70  mm,  fell  to  0.690  MPa  at  19.05  mm  and  remained  constant 
until  barrel  exit,  as  illustrated  in  Figure  7. 
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Figure  6.  Baseline  and  Chambrage  Configurations 


Hie  results  of  t^is  jecond  series  of  calculations  are  given  in  Table  6.  The  propellant 
dimensions  selected  for  these  calculations  were  those  determined  to  achieve  a  peak  pressure  of 
345  MPa  with  XKTC,  assuming  zero  barrel  resistance,  no  chambrage  and  the  propellant  ignited 
at  time  zero  for  the  given  C/M. 

Table  6.  Comparisons  of  Calculated  Maximum  Breech  Pressures 
and  Muzzle  Velocities  for  XKTC  and  IBHVG2 

Maximum 


Code 

Variation* 

Breech 

Pressure 

(MPa) 

Time 

(ms) 

Muzzle 

Velocity 

(m/s) 

Time 

(ms) 

C/M 

XKTC 

B 

345 

6.5 

1352 

10.4 

1.0 

IBHVG2 

B 

346 

6.3 

1363 

10.3 

1.0 

XKTC 

R 

363 

6.5 

1379 

10.2 

1.0 

IBHVG2 

R 

365 

6.3 

1387 

10.2 

1.0 

XKTC 

F 

332 

6.0 

1337 

9.8 

1.0 

XKTC 

C 

310 

6.5 

1303 

10.5 

1.0 

XKTC 

FRC 

321 

6.0 

1317 

9.8 

1.0 

XKTC 

B 

345 

11.2 

779 

18.1 

0.25 

IBHVG2 

B 

348 

11.1 

783 

18.0 

0.25 

XKTC 

R 

363 

11.2 

792 

17.9 

0.25 

IBHVG2 

R 

367 

11.0 

7% 

17.8 

0.25 

XKTC 

F 

347 

9.3 

779 

16.2 

0.25 

XKTC 

C 

336 

11.4 

774 

18.1 

0.25 

XKTC 

FRC 

351 

9.7 

784 

16.5 

0.25 

*  B  *  Baseline  R  =  Resistance  F  *  Flamespreading  C  =  Chambrage 

The  first  line  of  Table  6,  the  baseline  calculation  with  XKTC,  is  the  same  as  the  XKTC 
calculation  on  the  second  line  of  Table  5.  The  second  line  of  Table  6,  the  baseline  IBHVG2 
calculation,  is  also  the  same  as  that  shown  on  the  second  line  of  Table  5.  The  next  two  lines 
result  from  adding  bore  resistance,  and  show  the  expected  rise  in  peak  pressure  for  both  XKTC 
and  IBHVG2.  The  following  line  results  from  adding  just  flamespreading  to  the  XKTC  baseline 
calculation  and  documents  a  drop  in  peak  pressure.  The  next  line  represents  adding  just 
chambrage  to  the  baseline  XKTC  calculation  --  note  the  huge  drop  in  peak  pressure  of  35  MPa! 
The  final  calculation  in  the  C/M  =  1.0  series  has  flamespreading,  bore  resistance,  and 
chambrage  added  to  the  baseline  data  base.  Again,  the  large  decrease  in  peak  pressure  is 
attributed  primarily  to  the  chambrage. 

For  C/M  of  0.25,  it  is  seen  that  bore  resistance  makes  a  significant  change,  but  about 
the  same  for  both  XKTC  and  IBHVG2.  The  influence  of  flamespreading  and  chambrage  are 
about  2/5  as  great  as  those  for  the  higher  C/M. 
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Plots  of  ratios  of  breech  pressure  to  base  pressure  from  the  XKTC  calculations  for  this 
second  series  of  comparison  calculations  are  found  in  Figure  8.  Plot  25  results  from  having 
added  flamespread  alone  to  the  bassline.  Plot  26  has  resistance  added  to  the  baseline.  Neither 
Plot  25  nor  Plot  26  shows  much  difference  from  the  baseline  plot.  Plot  28  had  added 
chambrage  to  the  baseline;  Plot  29  had  added  chambrage,  flarnespreading,  and  resistance.  Both 
Plots  28  and  29  show  a  much  lower  pressure  ratio  at  early  time.  Plots  27  and  30  are  "full" 
XKTC  simulations  (which  now  include  projectile  intrusion  into  the  chamber)  of  the  inert  case 
(Plot  27)  and  the  caseless  gun  firings  (Plot  28).  The  fact  that  the  gradient  curves  in  Figure  8 
have  assumed  the  general  shape  of  the  measured  curves  leads  us  to  believe  that  the  major 
interior  ballistic  parameters  have  now  been  included. 


Figure  8.  XKTC-Calculated  Breech/Base  Pressure  Ratio  Curves  With  Added  Complexities. 
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V.  DISCUSSION  AND  CONCLUSIONS 


Earlier  studies  using  the  XKTC  code1  demonstrated  that  propellant  packaging  and 
boattad  intrusion  can  have  significant  impact  on  maximum  chamber  pressures.  The  current 
investigation  reveals  that  flamespreading  and,  to  an  even  greater  extent,  chambrage,  can  affect 
calculated  pressures  as  well,  particularly  at  large  C/M  ratios. 

In  this  study,  we  have  seen  that  the  presence  of  chambrage  makes  a  significant 
difference  in  the  pressure  gradient.  The  larger  cross  section  of  the  chamber  results  in  a  closer 
axial  proximity  of  combustion  product  gases  to  the  projectile  base,  with  less  axial  motion  and 
shorter  transit  times  required  to  transfer  pressures  downbore.  The  result  is  a  significant 
reduction  in  the  pressure  gradient,  as  shown  in  plots  28  and  29  of  Figure  8. 

The  role  of  flamespreading,  as  well,  has  been  demonstrated  in  this  study.  The  phased 
ignition  of  propellant  surfaces,  rather  than  the  simultaneous  ignition  event  assumed  in  most 
lumped-parameter  codes,  can  influence  the  overall  rate  of  gas  production,  as  well  as  the 
formation  of  pressure  waves,  resulting  in  differences  in  the  inbore  trajectory  and  impacting 
maximum  chamber  pressure.  Wave  dynamics  associated  with  the  rarefaction  accompanying 
projectile  motion  may  add  further  structure  to  the  pressure  gradient. 

But  a  more  interesting  feature  of  the  structure  of  the  pressure  gradient  is  shown  in  this 
study  to  accompany  the  formation  of  a  region  of  axial  ullage  between  the  propellant  charge  and 
the  projectile  as  it  first  moves  downbore.  XKTC  calculations  suggest  that  a  lower  gas  pressure 
gradient  exists  in  this  .gle-phase  region  than  in  the  two-phase  region  of  the  propellant  charge, 
primarily  because  of  a  lower  resistance  to  the  transfer  of  pressure  information.  This  result 
suggests  that  lumped-parameter  codes  might  benefit  from  a  two-region  pressure  gradient  in 
order  to  capture  the  true  structure  of  the  pressure  field.  Such  a  feature  might  prove  to  be 
particularly  important  for  simulation  of  stick  charges,  for  which  a  well  defined  boundary  between 
the  two  regions  persists  throughout  the  interior  ballistic  cycle,  or  for  highly  nonuniform  initial 
distributions  of  propellant,  as  when  firing  low-zone  artillery  changes. 

All  of  the  above  effects  are  exacerbated  by  an  increase  in  C/M,  since  the  projectile 
then  moves  out  more  rapidly,  and  axial  dime  ions  and  accompanying  transit  times  are  increased 
at  a  time  when  significant  amounts  of  gas  are  still  being  locally  produced  in  the  gun  chamber. 

We  conclude  from  this  study  that  lumped-parameter  interior  ballistic  codes  could 
benefit  greatly  from  the  inclusion  of  a  new  or  modified  gradient  equation  including  functional 
dependence  on  C/M,  chambrage,  propellant  distribution  and  ullage;  the  influences  of 
flamespreading  and  wave  dynamics  may  also  be  included,  though  the  basis  for  such  terms  would 
necessarily  be  more  heuristic. 
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APPENDIX  B 
IBHVG2  Data  Base 
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